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ABSTRACT 

This  is  a  report  of  a  search  for  propagation  anomalies  using  a 
large  quantity  of  high  frequency  data  produced  as  a  byproduct  of 
BRIGHAM,    a  Department  of  Defense  project.      The  BRIGHAM  data  is 
based  on  890  KHz  wide  samples   of  the  HF  spectrum  at  a  25  cycle  rate, 
using  a  2.  8  KHz  resolution  for  a  duration  of  approximately  2.  5  minutes. 
This  method  of  data  collection  is  unique  and  it  was  hoped  that  propaga- 
tion anomalies,    including  wide  band  anomalies,    might  be  detected. 
Anomalies  are  believed  to  occur  in  the  propagation  of  radio  signals 
and  they  are  usually  other  than  known,    routine  effects  but  may  include 
known  effects  which  cannot  adequately  be  explained.      The  scope  of  this 
examination  was  limited  to  the  visual  analysis  of  computer  processed 
data  presented  on  an  interactive  graphics  unit. 
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I.     OBJECTIVES 

This  project  provided  the  opportunity  to  examine  unique  data 
gathered  with  the  AN/FLR-15  wide-band  scanning  receiver  and  the 
AN/FRD-10  circularly  disposed  antenna  system.      The  project  was 
suggested  by  the  availability  of  this  large  volume  of  data  which  had 
been  gathered  at  great  expense  for  another  purpose.      The  data  seemed 
to  offer  great  potential  and  the  possible  basis  for  considerable  research 
and  understanding  about  electromagnetic  propagation.      The  initial  con- 
cept was  to  observe  data  in  the  different  form  of  a  3-dimensional  con- 
tour of  amplitude,    frequency  and  time.      This   concept  included  the  use 
of  a  computer  and  interactive  graphics  unit  to  vary  the  size  and  dim- 
ensions of  the  contour  to  obtain  a  new  view  of  the  high  frequency  region 
of  the  spectrum  and  perhaps  observe  propagation  characteristics  not 
easily  recognized  otherwise. 

The  objectives  of  this  project  were  to  take  the  given  sets  of 
data,    write  the  computer  programs  necessary  to  present  the  data 
visually  and  to  analyze  the  data.      The  objectives  of  the  analysis  were 
twofold:     To  study  propagation  effects  on  high  frequency  signals, 
becoming  familiar  with  the  data,    the  mode  of  visual  presentation 
and  the  appearance  of  routine  propagation  effects;  To  look  for  unusual 
and  anomalous  propagation  effects. 


II.     INTRODUCTION 

It  is  stressed  that  the  scope  of  this  effort  was   confined  to  a 
visual  analysis  of  the  data  for  detecting  unusual  patterns  and  anomalies 
which  could  be  attributed  to  propagation  conditions.      The  search 
was  conducted  in  the  high  frequency  region  (3-30  MHz)  using  data 
on  computer  cards.      The  data  cards  -were  produced  during  a  phase 
of  BRIGHAM,    a  project  conducted  for  the  Department  of  Defense  by 
Sanders  Associates,    a  private  corporation.     Each  data  set  represents 
the  output  of  a  swept-tuned  panoramic  receiver  and  preserves  amp- 
litude information  vs  frequency  and  time.      For  visual  presentation, 
the  data  were  processed  on  a  Scientific  Data  Systems  XDS  9300 
computer  and  displayed  on  an  Adage  AGT/10  graphics  unit  (figure   1). 
The  equipment  is  part  of  the  Electrical  Engineering  Computer  Center 
of  the  Naval  Postgraduate  School.     Interesting  sets  of  sweeps  have 
been  extracted  as  one  picture,    plotted  on  a  CALCOMP  563  plotter 
and  included  in  this   report.      Forty- six  data  sets  were  available  for 
examination,    each  consisting  of  7195  cards.     Several  data  sets  were 
examined  in  detail  with  the  intent  of  using  other  sets  for  verification 
of  observed  anomalies.      Recommendations  for  further  examination 
are  given  at  the  end  of  this   report. 
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A.  NATURE  OF  THE  DATA 

The  data  were  gathered  using  a  swept-tuned  receiver  and  either 
an  omnidirectional  or  a   12     beam  antenna.      Each  data  set  was  produced 
by  sampling  the  receiver  output  as  it  scanned  downward  through  an 
890  KHz  band.     Sets  were  gathered  at  different  base  (starting)  fre- 
quencies and  at  different  times  of  the  day.      The  sets   contain  the 
amplitude  (in  db)  of  the  318  contiguous  frequency  bins  (sampling  points) 
and  consist  of  7194  cards  to  represent  the  3597  scans.      The  initial 
card  in  each  set  contains  date,  time  and  parameter  information  (fre- 
quency,   etc.  )  for  that  set. 

B.  PRESENTATION  OF  THE  DATA 

Computer  programs  for  (1)  transferring  the  data  from  cards  to 
tape,    (2)  visual  presentation  of  the  data  and  (3)  extraction  of  the  data 
for  graph  plotting  are  included  at  the  end  of  this   report.      There  are 
two  computer  programs  for  presenting  the  data  and  they  reflect  the 
trade-off  between  the  width  of  the  displayed  spectrum  and  the  number 
of  scans  which  could  be  viewed  simultaneously.      This   restriction  was 
due  to  core  memory  size  in  the  XDS  9300  and  again  in  the  AGT/10. 

1 .  Anomaly  A 

This  program  was  used  for  viewing  the  data  as  a  contour 
and  required  appreciable  dimensions   in  time  and  frequency  to  be  of 
use.     It  was  initially  adapted  from  a  library  program  but,    along  with 


two  accompanying  subroutines,    was   changed  almost  entirely  in  adapt- 
ing it  to  the  BRIGHAM  data.     Anomaly  Plot  A  was  tised  with  this 
program  for  plotting  graphs.      The  data  were  displayed  as  shown  in 
figure  2  with  the  most  recent  scan  at  the  bottom  of  the  picture.      For 
each  picture  change,    a  new  scan  was  brought  on  at  the  bottom  and  the 
oldest  scan  moved  off  at  the  top.      The  scans  were  not  moved  up  in 
a  continuous  or  movie-like  fashion  but  all  shifted  upward  a  step  at  a 
time  to  replace  previous   scans.     Options  were  available  for  bringing 
on  more  than  one  scan  at  a  time.      The  amount  could  be  varied  so  as 
to  change  the  entire  picture  at  once  if  desired.     Scans   could  be  skipped, 
taking  every  fifth  or  tenth  scan  for  example,    effectively  decreasing 
the  receiver  sampling  rate  to  accentuate  long  term  or  slowly  changing 
trends  in  the  data.     Also,    the  time  axis  could  be  rotated  to  vary  the 
3-dimension  or  surface  contour  effect.     Other  features  are  described 
in  the  computer  program  section.     An  optimum  size  of  130  bins  and 
20  scans  was  chosen  for  this  program.     A  deficiency  in  this  presenta- 
tion was  that  amplitude  changes  with  time  were  not  easily  observed. 
To  overcome  this,    the  Anomaly  B  program  was  written.      The  scale 
at  the  bottom  of  figure  3  and  all  the  graphs  is    10  frequency  bins 
(28  KHz)  per  division.      The  leftmost  bin  in  each  group  is  labeled  with 
the  computer  variable  ISTRT,    for  the  integer  value  of  the   "starting" 
bin.      ISTRT  is  a  NAMELIST  variable  and  its  value  can  be  changed 
at  the  graphics  console.      It  enabled  the  viewer  to  look  at  different 
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portions   of  the  band.      For  identification,    bins  are  located  by  counting 
to  the  right  from  the  starting  bin  (frequency  decreasing  toward  the 
right).      The  vertical  scale  of  the  Anomaly  A  presentation  is  5  db  of 
amplitude  per  division  (referred  to  receiver  sensitivity).     Scans  are 
spaced  1/25   second  apart  in  time  or  multiples  thereof  when  scans  are 
skipped. 

2.  Anomaly  B 

This  program  was  used  to  examine  more  carefully  the 
changes  in  a.mplitude  among  a  number  of  signals  (see  figure  3).     At 
the  top  are  five  traces   of  signals   chosen  by  specifying  the  bin  number. 
The  traces  show  the  amplitude  history  of  each  signal  up  to  a  maximum 
strength  of  30  db  each,    beginning  with  the  present  base  scan  at  the 
left  and  extending  to  the  right,    showing  the  last  80  scans.      The  am- 
plitude was  displayed  after  taking  the  maximum  value  of  the  signal 
in  three  bins   (the  one  selected  and  one  on  either  side)  and  scaling 
it.     At  the  bottom  are  displayed  all  318  bins  to  aid  in  selecting  the 
five  signals  to  be  viewed..     Because  of  the  memory  size  trade-off, 
only  five  scans  are  shown  simultaneously.      The  vertical  scales 
differ  from  those  in  Anomaly  A.      Because  of  a  smaller  scale,    the 
divisions  are   10  db  apart.     Anomaly  Plot  B  was  used  to  plot  this 
display  on  paper. 
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III.     ANALYSIS 

A.         PRELIMINARY  OBSERVATIONS 

Before  beginning  the  visual  analysis,    it  was  useful  to  predict 
the  type  of  displays  the  data  could  yield.      These  preliminary  observa- 
tions helped  visualize  numerous  features  and  signal  variations  not 
due  to  propagation  so  that  they  could  be  recognized  and  not  be  inter- 
preted as  anomalies  to  be  checked  and  verified.      Examples  of  these 
features  included  the  following: 

1.  Transmitter  characteristics:     Turning  on  and  off;  different 
types  of  keying  which  affect  frequency  and/or  amplitude. 

2.  Receiver  characteristics:     Filters  for  notching  out  inter- 
fering signals. 

Other  preliminary  observations  predicted  limitations  on  the 
ability  to  observe  valid  propagation  characteristics  which  might  be 
present.      These  limitations  were  due  to  receiver  characteristics, 
the  nature  of  the  data  and  the  method  of  visual  presentation. 

1.  Limitations  of  the  receiver:     The  890  KHz  band  restricted 
the  ability  to  relate  the  effects  on  harmonics  in  the  megahertz  region 
and  the  ability  to  observe  variations  over  a  wide  band. 

2.  Limitations  due  to  the  data:     The  lack  of  resolution 
imposed  by  the  2.  8  KHz  wide  frequency  bins  was  expected  to  make 
it  impossible  to  measure  slight  variations   in  frequency.      Better 
resolution  could  provide  better  identification  of  the  types  of  signals 
present  for  the  purpose  of  selecting  ones  with  more  stable  transmitted 
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characteristics  to  permit  more  reliable  analysis  and  measurement. 
The  short  duration  (2.  5  minutes)  sample  provided  a  mountain  of  data 
but  was   expected  to  prevent  observation  of  some  propagation  charac- 
teristics which  require  a  longer  period  to  manifest  themselves. 

3.  Shortcomings  in  the  visual  mode  of  presentation  were  not 

foreseen  until  much  later  in  the  project.     At  first  only  a  general  con- 
cept of  the  type  of  presentation  existed  and  the  final  characteristics 
were  considerably  different.     In  fact,    Anomaly  B  was  developed 
only  after  much  analysis  and  frustration  with  Anomaly  A.      It  should 
also  be  noted  that  although  the  B  version  was  developed  last,    it  was 
not  a  replacement  and  both  programs  complement  each  other. 

B.  PROCEDURE 

After  reading  a  data  set  onto  tape,    the  data  were  analyzed 
using  the  Anomaly  A  presentation  to  examine  the  entire  length  and 
width  of  the  set.      The  object  at  this  stage  was  to  find  gross  anomalies 
in  the  3 -dimensional  contour  which  might  be  easily  recognized. 
Hypothetical  examples  hoped  for  might  have  been  any  of  the  following: 
A  rapid  change  in  one  region  only  of  the  contour;   similar  changes  to 
one  or  more  amplitude  humps   (signals  or  combinations  of  signals)  in 
separated  but  perhaps  related  areas  of  the  contour;  holes  or  the 
absence  of  signals  in  limited  bands  of  the  contour.      The  Anomaly  B 
program  was  then  used  for  more  critical  analysis  of  amplitudes  over 
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a  longer  period.      The  amplitude -time  traces  displayed  data  over  a 
period  four  times  as  long  as  in  Anomaly  A  (80  data  points  vs  the  20 
scans)  and  examined  it  from  a.  different  perspective  (abscissa  in 
dimensions  of  time  rather  than  frequency  to  examine  individual 
signals). 

C.  FINDINGS  AND  OBSERVATIONS 

Data  sets  4,    5,    8,    9,    15,    32,    34  and  147  were  examined  in 
detail.     Set  4  was  the  first  set  examined  and  received  considerable 
attention  while  being  \ised  to  develop  the  first  computer  program.     A 
routine  for  numbering  each  scan  and  later  displaying  the  number  of 
the  base  or  bottom  scan  on  the  screen  was  developed  and  became 
indispensable  in  helping  identify  the  location  of  anomalies   or  unusual 
data  for  stibsequent  observations. 

The  first  unusual  feature  noted  is   shown  in  figure  4  (upper  left 
hand  corner).      The  amplitude  hump  there  is  about  4-6  KHz  wide  and 
is   seen  to  move-downward  in  frequency  with  time  at  a  rate  of  about 
250  KHz  per  second.     An  effort  to  trace  it  through  other  signals  was 
unsuccessful.      The  left  hand  margin  represents  the  upper  limit  of 
the  band  (18  MHz),    therefore  tracing  it  at  higher  frequencies  was  not 
possible.      The  rate  of  frequency  change  was  too  fast  to  be  classified 
as  a   "whistler"  or   "dawn  chorus".      These  phenomena  are  only  found 
at  considerably  lower  frequencies  (several  kilohertz)  and  the  bandwidth 
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FIGURE  4.        Data  set  4 
ISTRT  1 
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Base  frequency  18  MHz 


14 


of  the  amplitude  hump  was  too  great  for  this  type  of  classification. 
This  was  the  only  occurrence  of  this  phenomenon  which  was  observed. 
It  may  have  been  a  transmitter  being  tuned  or  a  special  purpose  signal. 

A  feature  present  in  sets  4,    5,    8,    9,    15,    34  and   147  are  sudden 
discontinuities  across  the  entire  spectrum  of  the  set.      Examples  are 
shown  in  figures  5  through  7.      The  transitions  were  too  rapid  to  be 
attributed  to  propagation  effects  and  a  satisfactory  explanation  has 
not  been  found.     Late  in  the  analysis,    in  set  34,    several  pictures  from 
Anomaly  A  were  pieced  together  (figure  7).      The  figure  shows  an 
effective  shift  of  160  scans.     At  one  point  it  seemed  to  be  the  result 
of  a  missing  data  card.      Two  cards  are  required  to  describe  one  scan 
and  if  one  is  missing,    the  next  card  and  all  subsequent  cards  are  all 
shifted  back  160  scans.     It  would  give  the  appearance  on  the  left  half 
of  the  scan  being  exchanged  with  the  right  half,    as  it  appears  in 
figure  7.     A  second  missing  card  would  then  put  them  back  in  the 
right  order,    ending  the  "discontinuity".      Unfortunately,    this   simple 
explanation  was  not  correct  because  an  examination  of  the  data  in  the 
region  of  the  set  34  discontinuity  revealed  no  missing  or  out-of-place 
cards.      Examination  of  cards  in  other  sets  with  discontinuities  also 
revealed  no  missing  cards.      The  bad  data  therefore  must  have  been 
due  to  some  sort  of  receiver  or  recording  malfunction. 

Figure  8,    instead  of  showing  any  anomalies,    shows  a  set  of 
well  ordered,    continuous  signals  of  nearly  constant  amplitude  and 
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FIGURE  5.        Data  set  4 
ISTRT   1 
Base  scan  1423 
Base  frequency   18  MHz 
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FIGURE   6.       Data  set  4 
ISTRT   100 
Base  scan  949 
Base  frequency  18  MHz 
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FIGURE    8.        Data  set  4 
ISTRT   104 
Base  scan  30 
Base  frequency  18  MHz 
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shape.      This  is  contrasted  in  figure  9  by  an  irregular  contour.     In 
addition  to  variations  in  amplitude,    some  signals   seem  to  be  turned 
on  and  off  alternately.      The  region  between  bins   30  and  38  contains 
a  signal  regularly  occurring  only  once  in  nearly  every  three  scans. 
At  bin  60,    the  pattern  is  irregular  and  is  probably  an  on-off  keyed 
signal. 

When  figure   10  was  first  observed,    it  was  thought  that  a  genuine 
anomaly  has  finally  been  found.      Band  limited  noise  seemed  to  appear 
first  to  the  right  of  bin  160  and  then  spread  to  higher  frequencies 
before  subsiding.      The  presence  of  the  notch  seemed  connected  also, 
marking  the  lower  limit  during  the  first  scan  and  persisting  before, 
throughout  and  after  the  noise  in  time.     After  some  study  and  con- 
sideration it  was  realized  that  each  sweep  required  a  finite  amount 
of  time  and  that  the  onset  of  noise  appeared,    although  abruptly,    at 
midscan.     It  is  seen  to  grow  weaker  with  time  before  finally 
disappearing. 

Further  analysis  ■was  conducted  with  Anomaly  B  to  examine 
amplitude  variations   of  individual  amplitude  humps.      Examination 

was  begun  in  set  32  because  it  was  free  of  the  discontinuities  present 

o 

in  each  of  the  previous  sets.     Also,    set  32  was  taken  from  a   12     beam 

antenna  rather  than  the  omnidirectional  one  used  in  the  sets  from 

4  to  15.      The  narrow  beam  was  chosen  to  confine  the  signals  observed 

to  a  particular  region  of  origin.     Although  points  of  signal  origin, 
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FIGURE  9.       Data  set  8 
ISTRT   1 
Base  scan  60 
Base  frequency  5  MHz 
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FIGURE   10.      Data  set  8 
ISTRT  100 
Base  scan   1137 
Base  frequency  5  MHz 
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and  consequently  propagation  conditions,    would  still  vary  with  each 
signal,    the  variations  would  not  be  as   extreme  as  in  the  360     case. 

Figures   11,    12  and   13  were  extracted  from  set  32  and  compare 
the  amplitudes  of  signals  found  in  bins    19,    64,    101,    156  and   182. 
Each  of  the  three  figures  were  drawn  with  a  different  sampling  rate 
but  have  the  same  terminal  or  base  scan  (leftmost  points  at  the  top 
of  the  figures  and  the  bottom  scan  in  the  lower  section  of  the  figures). 
Figure   11  shows  every  other  scan  (the  variable  NSKIP  set  equal  to  1 
in  the  computer  program)  beginning  with  scan  2840,    covering  a  6.4 
second  period.      Figure   12  shows  every  tenth  scan  (NSKIP=10)  begin- 
ning with  scan  2200  (covers  32  seconds).   Figure   13  shows  every  35th 
scan  beginning  with  scan  200  (covers   112  seconds).      The  scans  at  the 
bottom  of  each  picture  represent  the  entire  data  spectrum  and  show 
the  last  five  scans  at  the  particular  skip  rate  noted. 

A  few  similarities   exist  between  the  traces  in  figures   11  through 
13.      Propagation  effects  had  been  expected  to  create  corresponding 
dips  in  amplitude  among  the  signals  of  similar  points  of  origin  (such 
as  those  shown)  but  they  were  expected  to  appear  over  longer  periods. 
On  the  other  hand  however,    the  rapidity  of  the  fluctuations  is  not 
unusual.     Ionospheric  surface  irregularities  and  layer  displacements 
could  account  for  them  as  described  in  reference   1.      David  and  Voge 
[Rcf.    2],    for  example,    report  that  horizontal  displacements  have  been 
measured  by  doppler  techniques  at  velocities  up  to  60-80  miles  per 
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hoiir.     At  the  top  right  in  figure   11,    traces  2  and  3  show  an  inverse 
relationship.      Traces  in  figure   12  indicate  that  the  higher  frequencies 
-were  affected  first  with  the  effects  shifting  lower  with  time.     It  is 
realized  that  since  the  signals  selected  for  observation  here  were 
not  equally    spaced  in  frequency,    the  straight  line  analysis  may  not 
be  particularly  valid.     Also,   because  of  the  small  number  of  similar- 
ities in  these  samples,    they  may  only  be  coincidental. 

Figures    14  through  17  were  also  taken  from  set  32  but  compare 
the  amplitudes  of  bins  35,    85,    167,    224  and  303.      The  sampling  rates 
are  shown  for  each  figure  as  well  as  the  terminal  scan  number  and 
duration  of  the  sample.      Figure   14  shows  no  anomalies  or  correlation 
among  signals  with  NSKIP=2.     In  figure   15,    a  periodic  relation  of  longer 
duration  seems  to  exist  between  traces   1  and  3.      The  propagation  con- 
ditions affecting  the  two  signals  are  similar  but  with  the  higher  fre- 
quency signal  (trace   1  at  the  top)  being  affected  earlier  in  time.      These 
effects  are  seen  again  in  figure   16  but  at  a  higher  skip  rate  (dips 
closer  together)  ■where  they  can  be  observed  over  a  longer  period. 
No  data  exists  to  verify  this  pattern  for  the  preceding  time  interval 
but  it  does  not  appear  in  the  following  interval  (figure   17  -  covers 
the  interval  from  scan  1600  to  number  3200  -  NSKIP=20  over  80  data 
points).     It  may  have  died  out  or  the  display  characteristics  may  have 
prevented  it  from  being  seen.      The  relationship  stands  out  because 
of  the  frequency  dependence  of  the  effects.      This  type  of  anomaly  is 
characteristic  of  what   was    c>:jm-<  Led  whin  the   project   was   begun. 
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Figure   18  shows  an  unusual  relationship  between  the  signals  of 
traces   1,    2  and  3.      Two  of  the  signals  experienced  some  sort  of 
attenuation  -while  the  signal  in  the  middle  was  enhanced.     A  similar 
inverse  relationship  was  observed  in  figure   11.      The  same  phenomenon 
as  in  figure  18  can  be  seen  in  figure   19  but  at  a  higher  skip  rate.     Also 
present  were  periodic  variations  in  amplitude,    particularly  in  trace  3. 
Periodic  variations  became  more  pronounced  in  trace  4  (figure  20) 
in  a  later  time  interval. 

In  addition  to  observing  those  groups  of  signals  already  noted, 
other  groups   of  signals  were  also  chosen  in  order  to  sample  a  wide 
range  of  possible  combinations.     Some  groups   spread  the  sample  over 
the  entire  data  spectrum  while  others  were  confined  to  a  narrow  range 
of  bins.     A  spread  of  low  and  high  skip  numbers  was  also  included. 
Specifically,    the  following  groups  of  signals  in  set  32  were  among 
those  examined  but  with  nothing  significant  observed:     Bins  3,    9,    16, 
21  and  25  with  NSKIP  equal  to  1,    4,    16  and  40;   Bins    135,    147,    158, 
168  and  177  with  NSKIP  equal  to  3,    12  and  45;   Bins   8,    13,    17,    22  and 
27  with  NSKIP  equal  to  2,    10  and  35. 
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IV.      CONCLUSIONS  AND  RECOMMENDATIONS 

Although  a  tremendous  amount  of  data  was  available  as  a  by- 
product of  project  BRIGHAM  and  a  large  quantity  of  it  was   examined  in 
detail  during  this  effort,    the  many  anomalies 'hoped  for  at  the  begin- 
ning were  not  seen.      Three  significant  phenomena  were  observed  and 
preserved  in  this   report.      The  first,    the  signal  decreasing  in  freq- 
uency with  time  was  almost  certainly  not  a  product  of  a  propagation 
anomaly.     Its  identification  might  prove  to  be  an  interesting  challenge, 
however.      The  second  phenomenon  showed  that  its  periodic  features 
were  frequency  dependent.     Although  periodic  fading  or  intensity 
fluctuations  were  noticed  at  various  times  while  examining  the  data, 
none  were  so  well  defined  nor  were  they  so  well  related  between 
signals  as  in  this  case.      Comparing  the  ionospheric  layers  to  a  cloud 
whose  features  are  constantly  changing,    it  was  not  surprising  to  see 
the  phenomenon  confined  to  a  short  period  of  time.      Likewise,    the 
irregular  features  of  the  ionosphere  would  allow  several  patterns 
to  occur  at  once,    superimposed  on  each  other  to  create  the  overall 
observed  effect.     Some  of  these  secondary  features  seem  to  show 
through  in  figure   16.      The  third  phenomenon  was  the  simultaneous 
enhancement  and  attenuation  of  signals.      The  seemingly  incompatible 
conditions  which  caused  this  may  have  only  been  due  to  multipath 
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effects.     It  is   strange,    however,    that  no  such  pronounced  effects 
■were  observed  before  or  after  the  anomaly,    particularly  in  trace  3, 
and  that  they  occurred  simultaneously  among  the  three  signals.      It  is 
expected  that  a  great  number  of  other  anomalies  exist  and  that  many 
of  them  are  preserved  in  the  BRIGHAM  data.     Approximately  one 
fifth  of  the  available  data  were  examined  and  due  to  time  limitations, 
only  one  fifth  of  that  amount  was  examined  with  Anomaly  B.     Analysis 
with  Anomaly  B  was  particularly  time  consuming  because  only  five 
signals  could  be  examined  at  once.     It  is  likely  that  other  anomalies 
are  observable  with  the  existing  programs  and  many  others  detectable 
with  modifications  and  different  analysis  tools.     Some  recommendations 
for  improving  the  analysis  techniques  are  as  follows: 

1.  Amend  the  Tape  Load  program  to  read  the  binary  number 
of  each  data  card  and  alert  the  operator  to  missing  cards.     A  one-up 
numbering  system  to  keep  pace  with  card  numbers  could  be  used 
with  an  "if"  statement  to  trigger  an  output  statement  when  a  card  is 
discovered  missing. 

2.  Modify  the  procedure  for  selecting  data  for  the  amplitude 
history  traces  in  Anomaly  B.     When  cards  are  skipped  to  examine 
long-term  propagation  effects,    an  averaging  routine  should  be  entered 
to  remove  the  short  term  amplitude  fluctuations.      The  presence  of 
the  short  term  effects  tends  to  mask  those  manifested  only  over  a 
longer  period  of  time.     A  modification  of  tlii:    type  was  attempted 
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near  the  end  of  the  analysis  but  was  unsuccessful  due  to  a  lack  of 
time  needed  to  perfect  it. 

3.  Incorporate  a  feature  for  backing  up  the  tape  by  a  variable 
(NAMELIST  variable)  number  of  scans.      Two  attempts  to  do  this 
failed.     In  one  case  it  was  due  to  the  fact  that  each  scan  constituted 
one  "physical  record"  in  length  and  that  backspacing  was  based  on 
"logical  records".     If  accomplished,    the  feature  would  be  a  significant 
time  saving  addition. 

4.  Temporarily  abandon  the  visual  analysis  approach  of 
Anomaly  B  and  employ  statistical  analysis  techniques.      Evaluation 
of  the  correlation  between  various  sets  of  signals  could  be  conducted 
rapidly  and  rather  exhaustively  on  a  larger  computer  such  as  the 
IBM  360  available  at  the  Naval  Postgraduate  School.      Preliminary 
analysis  with  Anomaly  A  might  prove  useful  in  the  selection  of  signals 
for  analysis  and  indispensable  in  locating  the  discontinuities  described 
earlier  so  as  to  avoid  misinterpreting  results.  /• 

5.  Use.  of  a  larger,    faster  computer  such  as  the  IBM  360  to 
load  cards  onto  magnetic  tape  would  significantly  reduce  "dead  time" 
in  an  extension  of  this  project.      Over  an  hour  was  required  to  load  a 
single  data  set  onto  tape  with  the  equipment  described  in  this   report. 
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GLOSSARY 

1.  Bin  (frequency  bin):     A  2.  8  KHz  wide  data  sampling  point. 
Energy  present  in  these  limits  was  qiiantized  and  represented  in  the 
range  from  0-28  db. 

2.  Data  set:     3597  successive  sampled  sweeps   or  scans   of  the 
receiver  and  preserved  on  punched  cards  of  magnetic  tape. 

3.  Dawn  chorus  or  whistler:     Interference  in  the  audio  range  from 
storm  discharges  or  high  energy  particle  interaction  in  the  ionosphere. 
Dawn  chorus  is  produced  around  sunrise  and  sounds  like  a  serenade 

of  birds.      It  is  a  rise  in  frequency  above  about  2  KHz.     Whistlers  are 
audible  decreases  in  frequency  and  are  propagated  from  storm  dis- 
charges along  magnetic  force  lines  of  the  earth  in  the  range  from  a 
few  hundred  hertz  to  20-30  KHz. 

4.  NAMELIST  variables:  Computer  variables,  properly  specified, 
which  can  be  entered  or  changed  after  the  program  has  been  compiled 
and  is  running.  These  variables  were  used  to  control  various  aspects 
of  the  program.      The  inost  used  were: 

ISTRT  -  Integer  number  of  the  leftmost  or  starting  bin  to  be 

displayed. 

IWDTH  -        Integer  number  of  adjacent  bins  to  be  displayed. 

NSKIP  _  Number  of  scans  to  be  skipped  over. 

NMBRJMP  -  Number  of  scans   to  be  brought  on  (and  moved  off) 
the  screen  at  one  time. 

5.  Scan:     One  sweep  of  the  receiver  from  one  end  of  the  band  to  the 
yther. 
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The  following  subroutines,    GINP,    GINPUT,    FNS  and  VCD,    are 
common  to  and  required  by  both  programs  Anomaly  A  and  Anomaly  B. 
They  are  library  subroutines  but  are  nonstandard  and  subject  to  change. 
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SUBROUTINE  "GINPUT" 
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SUBROUTINE  "FNS1 
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SUBROUTINE   "VCD1 
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PROGRAM  "GRAPH  PLOT" 

"GRAPH  PLOT"  is  a  library  program  on  paper  tape  and  is  used  with 
the  CDC  160  computer  to  drive  the  CALCOMP  563  plotter.     It  is 
found  in  the  Electrical  Engineering  Computer   Center  of  the  U.    S. 
Naval  Postgraduate  School. 
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